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Abstract. In this paper we present some new binomial identities for multiple harmonic 
sums whose indices are the sequences ({l} a , c, {l} 6 ), ({2}°, c, {2} 6 ) and prove a number of 
congruences for these sums modulo a prime p. The congruences obtained allow us to find nice 
p-analogues of Leshchiner's series for zeta values and to refine a result due to M. Hoffman and 
J. Zhao about the set of generators of the multiple harmonic sums of weight 7 and 9 modulo p. 
As a further application we provide a new proof of Zagier 's formula for C({2} a , 3, {2} b ) based 
on a finite identity for partial sums of the zeta-star series. 



1. Introduction 

In the last few years there has been a growing attention in the study of p-adic analogues 
of various binomial series related to multiple zeta values, which are nested generalizations of 
the classical Riemann zeta function £(s) = X)n=i V 72 *- The main reason of interest of such 
p-analogues is that they are related to divisibility properties of multiple harmonic sums which 
can be considered as elementary "bricks" for expressing complicated congruences. Before 
discussing this further we recall the precise definitions of such objects. 

For r G N, s = (si, S2, . . . , s r ) E (Z*) r , and a non-negative integer n, the alternating multiple 
harmonic sum is defined by 

it I \ rr s g n ( s i) fc ' 

H n (s 1 ,s 2 , . . . ,s r ) = 2^ 11 — 

l<fci<fc 2 <...<fc r <n i=l 

and the "odd" alternating multiple harmonic sum is given by 

sgn (si) ki 



0<ki<k 2 <...<k r <n i=l v ' 

If all s\,...,s r are positive, then H n (s\, . . . , s r ) and iT n (si, s 2 , ■ ■ ■ , s r ) are called multiple 
harmonic sum (MHS) and "odd" multiple harmonic sum, respectively. For s = (s\, . . . , s r ) G 
N r , it is also convenient to consider the multiple star harmonic sum (or non-strict MHS) 

S n {si, ■ ■ ■ ,S r ) = ^2 Til 

fvi * * * tvr 

l<fci<---<fc r <n 1 

The integers l(s) := r and w := |s| := YH=i \ s i\ are ca hed the length (or depth) and the 
weight of a multiple harmonic sum. By convention, we put H n (s) = H n (s) = if n < r, and 
H n {$) = H n {%) = S n {%) = 1. By {s\, s 2 , ■ ■ ■ , Sj} m we denote the sequence of length mj with 
m repetitions of (s\, s 2 , ■ ■ ■ , Sj). 

MHSs occur naturally in different areas of mathematics such as combinatorics, number 
theory, algebraic geometry, quantum groups and knot theory. For a long time, MHSs have 
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been of interest to physicists [2, 30], who realized these sums as Mellin transforms of some 
functions occurring in quantum field theories. The limit cases of MHSs give rise to multiple 
zeta values (MZVs): 

C(si, s 2 , ■ ■ ■ , s r ) = lim H n (s 1 ,s 2 , ■ s r ), 

n— ¥oo 

(*(s 1 ,s 2 , ... ,s r ) = lim S n (si,s 2 , ... ,s r ) 

n—¥oo 

defined for s±, . . . s r -i > 1 and s r > 2 to ensure convergence of the series. 

The earliest results on MZVs are due to Euler who elaborated a method to reduce double 
sums of small weight to certain rational linear combinations of products of single sums. In 
particular, he proved the simple but non-trivial relation 2) = £(3) and determined the 
explicit values of the zeta function at even integers: 

C ( 2m ) = HTjggZ? , )2m 
sv ' 2(2m)! v ' 

where £ Q are the Bernoulli numbers defined by the generating function 

oo 



k 

X 

< - i ^- ' >h ~V 

k=0 



It is easy to verify that Bq = 1, B\ = —1/2, B 2 = 1/6, B4 = —1/30, and B 2m+ i = for 
m > 1. The denominators of the Bernoulli numbers B 2m are completely characterized by the 
Clausen-von Staudt Theorem [20, p. 233], which states that 



for m € N, B 2m + - is an integer 



n 

p-l\2m 

Early results for the values modulo p of the multiple harmonic sums, when p is a prime greater 
than 3, go back to the classical Wolstenholme's Theorem [32]: 

flp-i(l) = (mod p 2 ), flp-i(2) = (mod p). 

Glaisher [8] in 1900, and Lehmer [21] in 1938, proved that even the multiple harmonic sums 
H p -i(m) modulo a higher power of a prime p > m + 3, are related to the Bernoulli numbers: 

, . = J P 2 Bp-m-2 (mod p 3 ) if m is odd, 

lip — ]_ [TTl J — K 

{^+iP B p-m-i (modp 2 ) if m is even. 

The systematic study of MZVs began in the early 1990s with the works of Hoffman [16] 
and Zagier [33]. The set of the MZVs has a rich algebraic structure given by the shuffle and 
the stuffle (harmonic shuffle or quasi-shuffle) relations. These follow from the representation 
of multiple zeta values in terms of iterated integrals and harmonic sums, respectively. There 
are many conjectures concerning multiple zeta values and despite some recent progress, lots 
of open questions still remain to be answered. 

Let 3w denote the Q- vector space spanned by the set of multiple zeta values C( s i> ■ ■ ■ ,s r ) 
with s r > 2 and the total weight w = si + - • -+s r , and let 3 denote the Q- vector space spanned 
by all multiple zeta values over Q. A conjecture of Zagier [33] states that the dimension of 
the Q- vector space 3™ is given by the Perrin numbers d w defined for w > 3 by the recurrence 

dw = d w - 2 + d w -z 

with the initial conditions do = 1, d\ = 0, d 2 = 1. The upper bound dim3wj < d w was proved 
independently by Goncharov [9] and Terasoma [28]. 
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It is easy to see that the Perrin number d w is equal to the number of multiple zeta values 
£(si,...,s r ) with si + • • • + s r = w and each Sj £ {2,3}. While investigating the deep 
algebraic structure of 3, Hoffman [17] conjectured that the MZVs C( s i> • • • > s r) of weight w 
with Sj £ {2, 3} span the Q-space 3«i- Very recently, this conjecture was proved using motivic 
ideas by Brown [4]. So the main problem which remains open is proving that the numbers 
((s\, . . . , s r ) with sj £ {2, 3} are linearly independent over Q. 

According to Zagier's conjecture, a basis for 3™ for 2 < w < 9 should be given as follows 
(see [31]): 

w = 2, d 2 = l, C(2); 

w = 3, d 3 = 1, C(3); 

w = 4, dt = l, C(2,2); 

w = 5, d 5 = 2, C(2,3),C(3,2); 

w = 6, d 6 = 2, C(2,2,2),C(3,3); 

w = 7, d 7 = 3, C(2,2,3),C(2,3,2),C(3,2,2); 

w = 8, d 8 = 4, C(2,2,2,2),C(2,3,3),C(3,2,3),C(3,3,2); 

w = 9, d 9 = 5, C(2, 2, 2, 3), C(2, 2, 3, 2), C(2, 3, 2, 2), C(3, 2, 2, 2), C(3, 3, 3). 

Brown [4] proved that part of this conjectural basis, namely the multiple zeta values of the 
form £({2} a , 3, {2} b ) can be expressed in terms of ordinary zeta values. The exact formulae 
were found and proved by Zagier [34]: 

K 



c«2y,3, { 2rt = 2 y>ir (( 2o 2 ; 2 ) - a - 2--)( 26 2 ; ,)) c(2r + i W * 



r) 



and 



A' 



C({2}«,3,{2} t 



- 2 E 



2r 
2a 



J r,a 



(1 



-2r^ 



2r 
26+1 



C(2r + l)H*(K-r), (1) 



where = a + 6 + 1 and 
H(n) = C({2}») = 



7T 



2» 



iT(n) = C*({2}") = 2(1 - 2 1 - 2 ")C(2n), n > 0. 



" (2n + l)!' 

The theory of multiple harmonic sums modulo p bears many similarities with the theory of 
multiple zeta values. In the early 2000s Zhao [35, 37] began to generalize the Wolstenholme 
theorem and Glaisher-Lehmer congruences to other multiple harmonic sums with special 
emphasis to the cases where the sums are divisible by a high power of a prime. Later, Hoffman 
[18] showed that the algebraic setup developed in order to study MZVs can be extended to 
deal with MHSs as well. In [19], Hoffman described the possible relations modulo p between 
MHSs of given weight not exceeding 9. 

For a positive integer w, let c w denote the minimal number of harmonic sums of weight w 
which are needed to generate all MHSs of weight w modulo p for p > w + 1. For the first few 
values of c w , Hoffman [19] obtained the following table: 



w 


1 


2 


3 


4 


5 


6 


7 


8 


9 










1 





1 


1 


2 


2 


2 



The value eg = 2 was calculated conditionally under the assumption of the congruence 



1 



>p-i 



1,1,1,6) = — B, 



54 ~ p - 3 ' 648 

which was first conjectured and verified for all primes 10 < p < 2000 by Zhao [35, Prop. 3.2]. 



+ 



1889 



B, 



p-9 



(mod p), 
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The purpose of the present paper is to study multiple harmonic sums of the form 

S n ({2r,c,{2} b ), S n ({l} a ,c,{l} b ), c = 1,2,3,.... 

We prove some new binomial identities for these sums and apply them to obtain congruences 
modulo a prime p for 

Vi({2}M,{2f), Vi({2}",3,{2f), Vi({l} a .2,{l} 6 ). (3) 

This allows us to give a new proof of Zagier's formula (1) for £*({2} a , 3, {2} b ) based on a 
finite identity for S n ({2} a , 3, {2} b ) and to formulate its finite p-analogue. In addition, those 
congruences of the finite sums (3) sharpen Hoffman's and Zhao's results on MHSs of small 
weights. Indeed, we show that in weight 7 the number of generators C7 equals 1, and in 
weight 9 we prove congruence (2), which implies the equality eg = 2 unconditionally. From 
[19, Section 10] and Corollaries 4.1 and 4.2 below, the set of generators modulo p for multiple 
harmonic sums S p -i(s) of weight w = \s\ < 9 when p > w + 1 in terms of Bernoulli numbers 
is as follows: 



w 


= 1, 


Cl 


= 0, 


0; 






w 


= 2, 


C2 


= 0, 


0; 






w 


= 3, 


C3 


= 1, 


S p - 


2) = B V . 


-3; 


w 


= 4, 


c 4 


= 0, 


0; 






w 


= 5, 


C5 


= 1, 


Sp- 


.i(l,4) = B V . 


■5; 


w 


= 6, 


C6 


= 1, 


S p - 


-i(l,l,4) = - 


.1 R2 . 

6 n p-3' 


w 


= 7, 


C7 


= 1, 


Sp- 


-i(l,6) =B V . 




w 


= 8, 


C8 


= 2, 


Sp- 


-i(l,l,6), S p . 


_i(l,4)5 p _i(l,2) =Bp^B p 


w 


= 9, 


C9 


= 2, 


Sp- 


-i(l,8) = Bp. 


_ 9 , s r i(i.2)Vi(M,4) = 



Hoffman [19] conjectured that all multiple harmonic sums S p -i(s) (or i? p _i(s)) can be written 
modulo p as sums of products of "height one sums" S , p _i(l, . . . , 1, 2h) (or H p -i(l, . . . , 1, 2/i)) 
with rational coefficients. Note that this conjecture has been confirmed only for sums of 
weight w with w < 9. Zhao [36] continued computations of sets of generators for w = 10, 11, 12. 
Using reduction, stuffle and duality relations he found that the MHSs of weight 10 can be 
generated by the set 

Bp-zBp-i = Sp-\{\,2)Sp-i{\,6), Bp_ 5 = S'p_ 1 (l,4), 

S p _i(l,l,8), S p _i(l,l,l,l,6), Vi(2,2,l,4,l), 
which according to Hoffman's conjecture still contains the extra term ^,_i(2,2,l,4,l). Zhao 
[38] kindly communicated us that using our Theorems 4.1 and 4.2 (see below) it is possible 
to reduce the set of 8 generators for MHSs of weight 11 found in [36] to the set of 5 elements: 

Bp-n = S p -i(l, 10), B p - 5 Bp_ 3 , B p -3Sp-i(l, 1, 6), S'p_i(l, 1, 1, 8), ^-1(5, 3, 2, 1). 

In contrast to the corresponding conjecture for MZVs, not all "height one sums" are claimed to 
be linearly independent modulo p, because they satisfy the duality relation [19, Theorem 5.2]: 

Sp- 1 ({l} k - 1 ,h) = (-l) k+h Sp- 1 ({l} h -\k) (modp), p>max(M)- 

In [36, Remark 2.3], Zhao conjectured that weight-8 generators S l p _i(l,l,6) and B p ^B p ^ 
should be linearly independent modulo p over Q(10), where 

Q(w) := {a/b G Q : a/b is reduced and if a prime q\b then q < w}. 

Zhao [36] also formulated a general conjecture on linear independence modulo p of products 
of Bernoulli numbers 

Bp—i 1 Bp^i 2 ■ ■ ■ Bp—i r 
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over Q(w + 2), where odd indices greater than 1 with i\ + . . . + i r = w, which 

in particular, implies the linear independence of weight-9 generators B^_ 3 and -B p -g modulo 
p over Q(ll). 

As an application of our results to ordinary zeta values, we consider generalizations of 
Apery's famous series 



used in his irrationality proof of £(2) and C(3)- Among different extensions [22, 3, 1, 11, 13] 
of series (4), we focus on those obtained by Leshchiner [22] for values of the Riemann zeta 
function ((s) and Dirichlet beta function (3(s) = Y^ =Q {-l) n / {2n + Note that the value 
(3{2m + 1) as well as C{2m + 2) can be expressed in terms of tt, 

B(2m + D- (~ 1 ) m ^ 2m 7r 2m+i 
PVm + i)- 22m +2( 2m )! 7r 

where m is a non-negative integer and E2 m £ Q are Euler numbers. 

Leshchiner established the following four expansions in terms of multiple harmonic sums 
that generalize series (4): 

oo 



k — 1 V k ) 



(5) 



3=1 *=i ^ +2 ( 2 ^) 



■ C (2m + 3) = (-D^-iCm") 

{-l) k ~i H k ^({2} m -i) 



(6) 



1/ 2 2 m+2y K^™ + 2) " 4 JL l 6 fc(2fc + l) 2 

m oo^ (_l)fc+^i(2fc)ff fc({2} m-r 

+ Z^1> 16 fe (2fe + l) 2 J+ 2 

3=1 fc=o v ; 



(7) 



M7m ^ 3 - (-ir( 2 fc fc )g fc ({2p ^^ (-lr-^-fflgfc^r-j) 

/3(2m + i ) __2_ > 16fc(2fc + 1) + 1.1. I6*(2fc + 1)^+1 • (8) 

k — j — 1 k — 

Indeed, if we put m = in (5), (6), we get Apery's series (4). 

Recently, the authors [26, 14, 15], showed that the series (5), (6) for m = 0, 1 and the series 
(7), (8) for m = admit very nice p-analogues. Indeed, if we truncate the series (5), (6) with 
m = or 1 up to p — 1 and the series (7), (8) with m = up to (p — 3)/2, where p is an 
odd prime, we get congruences for the finite sums modulo powers of p expressible in terms of 
Bernoulli numbers. In this paper, we extend these results to any non-negative integer m. This 
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is done by showing that the finite p-analogues of the series (5)-(8) are related with MHSs of 
the form 

flp_i({2r,3), fl p _ 1 ({2r,l), ^E_i({2} m ,3), H P _i({2} m ,l). 

2. Identities for multiple harmonic sums 

We start with a list of binomial identities that we will need later. Note that a generalization 
of (13) appears in [10, (4.20)] 



Lemma 2.1. For any positive integers m,n and a non-negative integer I we have 



± (-d-'Ct*)-'-" 

k=i+i v y 

n i„(n\ „ In— 1 



l ( mn — 1 
n-l-1 



(9) 



Ifl>l, then 



If n > 2, then 



2T M2- = <ri (10 ) 

Z_> (n+k\ ln+l\ ' y } 

k=l+l \ k I \ I I 

n (n\ n -. 

^Ej% = Ei (id 

k=i K \ k ) k=i 



n lk\ ln\ 

V {l> = {l> (12) 

kkH k r) n n rr 1 J 

k=l \ k ) 



Proof. For the proof of our identities it is useful to consider the usual binomial coefficient (£) 
in a more general setting, that is, to allow an arbitrary integer to appear in the lower index 
of (£) . For this purpose we set 

k) [0 if integer k < 0. 

To prove the first identity we first observe that 



4 ) ( n _fe) =G(n,k + l)-G(n,k) (14) 



for positive integers n, k with 

G(n,*) = n _ fc 
Then summing both sides of equation (14) over k from I + 1 to n we obtain 



£ ("I)*- 1 ( ™ ) = G(n, n + 1) - G(n, / + 1) = -G(n, 1 + 1) = (-1)' ( ™ _M . 

I 7_1_1 \ / \ / 



k=l+l 

Similarly, for the proof of the second identity we have 



2k ( n ) 

-^. = G(n,k + l)-G(n,k) (15) 



(n+k\ 
\ k ) 
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for positive integers n, k with 

Summing both sides of equation (15) over k from / + 1 to n we easily obtain the result. 
To prove identity (11) we set 

f(m,k):= } k > and G(m,k) := 1 k \ -r- , m > 0, k > 1. 

fc( m + fe ) V ' 2(m+l)( m + fe ) 

Then it is easy to see that 

/(m + 1, k) - f(m, k) = G(m, k + 1) - G(m, fc). 

Summing both sides of equation (16) over m from to n — 1 we have 

n-l 

/(n,fe) = ^(G(m,Hl)-G(m,i)). 

m=0 

Now summing the above equation once again over k from 1 to n we obtain 

n n—1 n n—1 

E /(", k) = J2 E( G ( m ' k + X ) - = E ( G ( m ' « + 1) - G(m, 1)) 



fc=l m=0 fc=l m=0 

n—1 n—1 



^ v ' 1 2 ^ m + 1' 

m=0 m=0 



as required. 

For proving identity (12), it is easy to see that 

JL^ = G(k + l,l)-G(k,l) 

for positive integers k, I with 

Then summing both sides of equation (17) over k from I to n we get 

12 E ^TIJ = G(n + 1, Z) - G(Z, = G(n + 1, 1) = ^ 

To prove (13) we note that for integers n > 2, k > 0, 

,n+* =G(n,fc + l)-G(n,fc), 
I fc J 

with 

G(n - fc> - 2 ( »-i)m ' 

Then summing both sides of equation (18) over k from 1 to n we obtain 



E =G(n,n + l)-G(n,l) = 0, 

fc=i v fc / 



and the lemma is proved. 
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Lemma 2.2. Let k,m,n be positive integers, = (— l) k (™™ k )c^ where is an 

trary sequence independent of k, and a be a non-negative integer. Then for each c € N, 

1 " g fc -l(b);4ff _ " fffc-lCb)^ | ^ " g fc - 1 (b,B)^g 

n c /.a Z> £a+c m fcj 

fc=l fc=l j + \s\=a+c k=l 

j>0,s\>a 

where s = (s\,S2, - ■ ■ , s r ) G N r /or r > 0, and |s| = YH=i s i> K s ) = r - 
Proof. We first note that by (9), 

,(m) 



m (m) _ / 1 _ l\ .(m 



We prove (19) by induction on c. By (20), we have that 



n 



k=l 1=1 k=l+l 

= £H l - 1 (b)A%> (X i 



l a \n I 

i=i v 

1 " H^MAM " ^_x(b)Aff 

n 2^1 \a 2^1 la+1 

1=1 1=1 

which is (19) for c = 1. Now let c > 1. Then by induction, we have 
1 " g fc _ l( b)4ff 1 /l- ff^bj^f 

n c Z^ fca n c-l \ n 2^1 fra 

k=l \ k=l 

fc=i fc=i 



Z-^ fc a + c Z^ Z-^ 

fc=l " j + |s[=a+c fc=l 

j>0,si>a+l 

" H^a + 1)A^ ^ g fc - 1 (b,a + 1,8)^ 

+ m > : — = ; h m > m l(S> > — — 

k=l j"+|s|=c-l k=l 

j>0,si>0 

_ " ^_ l(b) Ag " g fc -i(b,s)<g 

Z^ /e a + c Z-/ Z^ 

fe=l j + |s|=a+c fe=l 

j>0,si>a+l 

" g^,(b,.)<g 

jr'+|s|=a+c k=l 
j>0,si=a+l 
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Theorem 2.1. For a positive integer n and non-negative integers a, b, and c > 2, 



^_ 1 (2a + i,s)(-l) fc - 1 (") ( 21 ) 



+ 4 V 2^ V 

Z^ Z^ i,2b+j(n+k\ 

i+j+\s\=c k=l \ k / 

«>lj>2,|s|>0 

where s = (s±, S2, ■ ■ ■ , s r ) £ W for r > 0, and |s| = Yli=i s «» ^( s ) = r - 

Proof. The proof is by induction on n. For n = 1 we have Si ({2}°, c, {2} b ) = 1, and the 
formula is true. For n > 1 we proceed as follows. If c = 2 then we should prove that 

^« 2 n = 2Eypf (22) 

/c — 1 V k ) 

where m = a + b + 1. First note that 

m 1 

l^ki^fe^-^k^n " lft2 Km /=0 n 

Then by the induction hypothesis, we have that 

n 1 n ~ l ( —-\\k-l (n-l\ 9 n-1 /_i\fc-l/'n-l\ m 2 Z 

c f/o\"i\ _ 9 1 \ " I L > \_k_l _ f \ ^ I I fc J \ ^ ^_ 

J ' Z_/ „2(m-0 Z_/ ^i^+fc-lA n 2m /n+fc-l\ 2-^^21 

1=0 k=l ^ V fc J fc=l v / Z=0 

2 ^ (— (£) (n 2m + 2 - fc 2m + 2 ) 



n 2m+2 Z^ £,2ra(™+^ 
ZZ- fc2m( n + k \ n 2m+2 ZZ/ /n+fcN 



and formula (22) easily follows by the equation (13). 
To prove (21) for c > 2 we note that for n > 1, 



1 



S n ({2}V,{2} b ) = E fc^W fc 2 -..fc 2 
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Then by the induction hypothesis and formula (22), we have 

b b 1 n-l (_l)fc-l("-l) 

S n {{2} a , c, {2} b ) = 2 , 2(a+0+c /„+fc-i\ 

Z=0 fc=l ft V k ) 



n—1 



i7 fe _ 1 (2a + z,s)(-l) fe - 1 ("- 1 ) 



+ 4 V - V 2^V 

^ Z^ „ 2(6-0 ^ ^ >2i+7 

«=0 i+j+ \ s \ =c k=l \ k ) 

i>l,J>2,|s|>0 

2 " (_i)fc-i^ 



k 



n 2b+c Z_> ^2a ( n + fc ) 



Changing the order of summation and summing the inner sum 

(V) £ « 21 _ ™ 2fe+2 - ^ 2b+2 (V) _ f« 26 ^ 2 



(«-i+ fc )Z^fc 2Z (n 2 -k 2 )k 2b ( n ~ l k +k ) \k 2b nV( n + fc )' 
we obtain 



n (_i)k-l(n\ 

S n ({2rM2} b ) = 2 ^J b+c{ ^ 



+ 4 V 2 , ( , 1 A g t - 1 (2a + i,s)(-l)*- 1 (;) 

2-~i (n+k\ 
i+j + \s\=c k=l \ k ) 

»>lj>2,|s|>0 

, _J_ v (- 1 ) fc ' i a) _ _L_ v (*) 

" " Z_> t,2a ( n + k \ n 26+2 Z^ i,2cH 



n 26+c 

k- 



-' fc2a( n + k \ n 2b+2 tfla+c-2 fn+k\ 

-1 \ k J k'=.\ \ k J 



4 V 2 *(s) V ^-i(2a + i,s)(-l) fc - i q) 

„26+2 Z-> Z-> hi-2(ri+k\ 

i+j+\s\=c k=l K \ k ) 

«>lj>2,|s|>0 

The final result follows as soon as we show that 

n c-2 Zw ^2a Zw £.2a+c-2 Z-^ Z_^ y^j-2 

fc=l fe=l j+j + | s |=c fe=l 

i>lj>2,|s|>0 

where A^J. = (-l) fe_1 (£)/("£ fe ) ■ This identity holds by Lemma 2.2 because 

2 £ 2t f)^ ^-i( 2 ° + M)<l _ £ oiM^ ^-iW^i , 

i+j+|s|=c fc=l jr'+|s|=2a+c-2 fe=l 

«>lj>2,|s|>0 i>0,si>2a 



□ 
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Corollary 2.1. For a positive integer n and non-negative integers a,b, 

1 I)*- 1 © 



1.2a (n+k\ 
k=l K \ k ) 



, il ^^*-i(2a + l)(-l) fc - 1 (2) 



(23) 



3„^J- , <5, 1Z| ) - I ^ fe2 (a+6)+3 (n+kj + 4 Z_> fc 2 b+ 2 (n+k) 

Theorem 2.2. For a positive integer n and non-negative integers a, b, c, 

wr,,M-±^+ E ± B ^' + ™-™ <"> 

fc=l j+j+|s|=c fc=l 

i>l,J>l,|s|>0 

Proof. Note that the case c = 1 is well known (see [6, Corollary 3]). It also follows from the 
polynomial identity 

E(l + x) kl - 1 _ ^ fn\ x k 
fei ■ ■ ■ fe m ~ 2-*>\k)k™ 

l<k 1 <...<k m <n 1 m k=l V 7 

which appeared in [27, Lemma 5.5]. For c > 1 the identity can be proved by induction on n. 
By a similar argument as in the proof of Theorem 2.1, we have 

s n ({ir, c , {i} b ) = J2 ^rs B _i({i} a , c , {iy) + ^s n ({m 

6+c / j fca / / f^f 1 — I ' J fca+c+l 



k=l 1=0 k=l 



,f_L V ^ g fc _ 1 (a + z,s)(-l) fc - 1 (V) 

Z^ n 6-i Z^ Z^ /ji+j 

Z=0 i+j+|s|=c fc=l 

i>lj>l,|s|>0 



Changing the order of summation and summing the inner sum 

6 

E 



1 \ & / / ?> 7 \ / 

n — 1 \ v-^ rr / n k\n 



k J ^— ' k l \ k b n \k 

1=0 



we obtain 



s,m^--±^ + E E gt - (B+ y c) 

fc=l j+j+|s|=c fe=l 

*>lj>l,|s|>0 

k=l ' k=l 

J_ v " ^(q + i.sX-l)*- 1 ® 

n 6+i z_> z^ jfcj-i 

«+j + |s|=c fc=l 
i>lj>l,|s|>0 
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Now the result easily follows from Lemma 2.2 with A^ k = ( — l) k ~ l (^ and the equality 

v A g ifc -i(a + i,B)(-l)*- 1 (;) v " jfr-xfrX-l)*- 1 ® 

Z> Z_> fc?'- 1 Z^ Z> # 

j+j+|s|=c fc=l j'+|s|=a+c-l fe=l 

»>lj>l,|s|>0 j>0,si>a 

□ 

Theorem 2.3. Lei a, 6 6e integers satisfying a > 1, 6 > 0. T/ien /or any positive integer n, 
S„(l,{2}') = 2g ;i J| s , (25) 

*.({2}M, { 2>>) = 2 vj JiX'ft - 4 E g ^."ff ' < 26 > 

/c — 1 V k } — 1 V k ) 

Proof. To prove the first identity we proceed by induction on b. For b = its validity follows 
from (11). Assume the formula holds for b > 0. Then we easily obtain from (12), 

' (?) 



£2 ff.2 ]2m+l(k+l) 

k=l k=l 1=1 \ I I 

n ^ n ik\ n (n\ 

\l) _ o \l) 



2 'Y2 pm+l k 2 (k+lj 2 pm+3 (n+lj ' 

as required. To prove the second identity of our theorem we proceed by induction on n. 
Obviously, it is valid for n = 1. For n > 1 we use the equality 

S n ({2r, 1, {2} b ) = £ -A-^ n _ l( {2}», 1, {2}') + -i I 5 n ({2}») 
z=o n 

and apply the same arguments as in the proof of Theorem 2.1. Then with the help of the 
formula (10) we may easily deduce the result. □ 

Remark 2.3.1. Note that the identities for the finite multiple harmonic sums found in The- 
orems 2.1, 2.3 can be used for evaluating multiple zeta-star values. For example, letting n 
tend to infinity in (23), (25), (26) we get expressions for C*({2} a , 3, {2} 6 ), C*({2} a , 1, {2} b ) in 
terms of alternating Euler sums: 

C({2} a , 3, {2} h ) = 2((2a + 2b + 3) + 4 £ t}^^^l±A , (27 ) 

k=i 



where 



C({2} a , 1, {2} b ) = 2C(2a + 2b + 1) - 4 ]T gfc , 2 \ ( +1 2q) , a,6>l, (28) 

fe=i 

C*(1,{2} 6 ) = 2C(26+1), 6>1, (29) 
00 / -i \fe— l 

cw :=E L ^- = ( 1 - 2l "*)cw. 



fc=i 



TOi/t £(1) = log 2, is t/ie alternating zeta function. Formula (29) was proved previously by 
several authors (see [23, p. 292, Ex. b], [39], [29]). Evaluations for length 2 Euler sums of 
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odd weight appearing on the right-hand sides of (27) and (28) in terms of zeta values are well 
known (see, for example, [7, Theorem 7.2]). Therefore by [7, Theorem 7.2], we obtain 

<W,M2>») = *E(( 26 + ,) - +v. - (£))«* + * <2A " " 2r) - 



where C(0) := 1, K = a + b + 1, which gives another proof of Zagier's formula (1) for the 
zeta-star value C*({2} a , 3, {2} b ), and 

C*({2}M,{2} b ) = 4f; 



r=l 



2r 



2a 



l) i 1 " i)) C(2r + 1)C(2a + 26 " 20 (30) 



/or a, 6 > 1. Note that formula (30) was also proved in [25, Theorem 1.6] by another method. 

3. Auxiliary results on congruences 

In this section we collect several congruences that will be required later in this paper, 
(i) ([35, Theorem 1.6]) for positive integers a,r and for any prime p > ar + 2, 

'(-l) r SS§P 2 £ P -ar-2 (modp 3 ) if ar is odd, 
{-^) r ~ l -^+lPB p -ar-i (mod p 2 ) if ar is even; 



flp-l(W) EE 



(ii) ([35, Theorems 3.1, 3.2]) for positive integers a±, a 2 and for any prime p > a\ + a 2 , 

(_1)«2 ^ ai + a2 > 

Ol 



Hp-i (ai,a 2 ) = 



B 



p—ai—a2 



(mod p), 



ai + a 2 

moreover, if ai + a 2 is even, then for any prime p > ai + a 2 + 1, 



flp_i(ai,a 2 ) e= p 



(-l) ai a 2 

Bp— ai -a 2 — l 



ai + a 2 + 1 
ai 

(mod p 2 ); 



-(-l) ai ai 



ai + a 2 + 1 

02 



— ai — a 2 



2(ai+a 2 + l) 

(hi) ([35, Theorem 3.5] if a±, a 2 , 03 G N and w := a\ + a 2 + 03 is odd, then for any prime 



p > w, we have 

#p_i(ai,a 2 ,a 3 ) = 



(-1) 



ai 



(-1) 



«3 



B 



p—w 



2w 



(mod p); 



(iv) ([24, Theorem 5.2]) for a positive integer a and for any prime p > a + 2, we have 
'-2g p (2)+^ (2 )_ p 2(2 9 3 (2) + ^ B?) _ 3 ) (modp 3 } ifo = 1> 

if a is even, 
if a > 1 is odd, 

where q p (2) = (2 P_1 — l)/p is the so-called Fermat quotient; 



Ht-i(a) 



^+TT p ( mod P 2 ) 

_2^2 Bp _ a ( m od p) 



(v) ([15, Lemma 1]) if a, b are positive integers and a + b is odd, then for any prime 
p > a + 6, 
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(vi) by (i) and (iv), for any positive integer a and any prime p > a + 2, we have 
flp-i(-a) = 2 1 - a H E _ 1 (a) - F p _i(a) 

'-2a p (2)+p< 7 2 (2)-p 2 (fg3(2) + \B p ^) (mod p 3 ) if a = 1, 



a ^ 2 — -pB p _ a -i (mod p 2 ) if a is even, 



o+l 

, ~ 2< ' 1 ~a — ~ B p - a (mod p) if a > 1 is odd; 

(vii) ([27, Theorem 3.1]) for positive integers a, 6 of distinct parity and a prime p > a+b+1, 

iJ p _i(-o, b) = Hp-^a, -b) = — — — — B p ^ h (mod p), 

g p -i(-o,-6)= Q + b (-l) b ( | K_ B -*> (modp); 

(viii) ([27, Theorem 4.1]) if a, 6, c £ N and iw := a + 6 + c is odd, then for any prime p > w, 

2H p -i(a, -b, -c) = H p -i(c + 6, a) + i7 p _i(— c, -6 - a) 

- Hp-i(-c)H p -i(-b, a) (mod p), 

2H p - 1 (-a,b, -c) = i/p_i(-c)iJ p _i(6, -a) - iJ p _i(-c, b)H p _i(-a) 

+ H p -i(—c — 6, —a) + iJ p _i(— c, —6 — a) (mod p). 

In [5, Lemma 6.2] it was proved that for positive integers a, 6 and a prime p > 2o + 26 + 1, 

flp_i(-2a, -26) = (mod p). 
However, we will need a stronger version of this congruence. 

Lemma 3.1. Zei a, 6 6e positive integers and a prime p > 2a + 26 + 1. T/ien f/ie following 
congruence holds modulo p 2 : 

, „ /(a-6)(l-2- 2a - 2b ) /2a + 26\ a + 6 \ 

flp-i(-2a, -26)=^ (2a + 1)(2fc+1) ^ 2fl J-2aT26TlJ^- 2a - 26 - 1 - 

Proof. Clearly (see [27, Lemma 7.1]), 

(-iy+ k „ (_i)p-j+p-fe 

H p -l(-2a, -26) = 2^ ~~2afJ2b~ = („ _ 7) 20(2, _ M26 

i<j<fe< P -i J i<fe<j<p-i VF 1/7 v ^ 7 

i<fc<i<p-i 

+ 2bpH p ^i(-2b - 1, -2a) + 2apH p _ 1 (-2b, -2a - 1) (mod p 2 ). 
Therefore by (vii), we obtain 

flp_i(-2a, -26) = flp_i(-26, -2a) 

2p(a-6)(l-2~ 2 ^ 26 ) /2a + 26\ , 2 (31) 

+ (2a + l)(26 + l) - ( 2a ( mod ^ 
On the other hand, we have 

flp_i(-2a, -26) + flp_i(-26, -2a) 

= flp_i(-2o)flp_i(-26) - flp_i(2a + 26) (mod p 2 ) 
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and therefore by (vi) and (i), 

flp_i(-2a, -26) + flp_i(-26, -2a) = -# p _i(2a + 26) 

(2« + 26)p 2 (32) 

'" 2a + 26 + l ^- 2a - 2b - 1 (m ° dp) - 
Now from (31) and (32) we obtain the required congruence. □ 

Lemma 3.2. For positive integers a, 6 and a prime p > 2a + 26 + 1, £/ie following congruence 
holds modulo p : 

2 2a+26_ 1 / x /2a + 26+l\ \ 

H P -i -2a, -26 - 1) = ; „ -rj— +1 B _ 2o -26-i- 

V v ' ; 2a + 26 + 1 I 2 2a+2b+1 V 2a / / p 



Proof. It is easily seen that 

fl p _ 1 (-2a, -26-1) -fT^-i (-2a, -26-1)= £ 1 ^ ^ 

(P-l)/2 / i No-fc rr /ox (P-l)/2 



(33) 



(-ir fc # p - fc -i(-2a) lP " (-l) fc ff p _ fc _ 1 (-2a) 
= L (p-fc^H-l 55 Ij ( mod ^- 

For the harmonic number # p _fc_i(— 2a), by (vi) we have 

H p -k-i(-2a) = ^2 = ^2 - = - ^2 - 2a 

j=l J j=k+l KF J > j=k+l J 

= -(flp_i(-2o) - H k {-2a)) = H k {-2a) (mod p). 
Substituting the above congruence in (33) we obtain 

(p-l)/2 ( _ ukf r ,_ s 

fl p _ 1 (-2o,-26-l)-^(-2a,-26-l)= 5] 1 ^r£S 

fc=i 

= #£^I (-2a, -26 - 1) + #2-1 (2a + 26 + 1) (mod p) , 

2 2 

which by (iv) and (vii) implies the required result. □ 

4. New congruences for multiple harmonic sums 

In this section we state and prove our main results on multiple harmonic sums. 

Theorem 4.1. Let a, b be non-negative integers and p be a prime such that p > 2a + 26 + 3. 
Then 

it o r 2lh _ (-l) a+6 (a-6) /2a + 26 + 2\ 

# p _!({2} ,3,{2}) = (fl + 1)(6+1) ^ 2a + 1 )B p -2a-2 b -S (modp). 
Proof. From Theorem 2.1 we have 

c fi2ni2lV2 r (-^(V) ^ g^a + DC-i)*- 1 ^ 1 ) 

^p-lU^ ,d,^J. j - fc2a+2b+3 / p -l+fcN + 4 Z : > A .26+2(p-l+fe) 

/c — 1 V k ) k — 1 V k ) 
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Since 



(V) _ (p-i)(p-2)...(p-fc) _ jj (i _ ^ fi(i+ p 



(f-i+ fc ) p(p + 1) ■ ■ ■ (p + A; - 1) p Al ^ j ; 11 ^ j 

ee <=^*(1 -pff*(l))(l -pff*-i(l)) ( 34 ) 



(-l) fe A; 



we obtain 

p-i 



(l-2pfT fc _i(l)-|) (modp), 



Vi({2} B , 3, {2} b ) S £ (l - 2p^_ l( l) - |) 

^ k=l 

■ ^ ff fc _ 1 (2g + l) / p\ 
• =i 

H p -\{2a + 26 + 2) + 4flp_i(l, 2a + 26 + 2) 



^ fc=i 



p 

- -H p ^(2a + 1, 26 + 1) + 8H p ^(2a + 1, 1, 26 + 1) 
p 

+ 8flp_i(l, 2a + 1, 26 + 1) + 8# p _i(2a + 2, 26 + 1) 
+ 4flp_i(2a + l,26 + 2) (mod p). 

Now by (i)-(iii), we easily obtain the required congruence for S p -i({2} a , 3, {2} b ). To prove 
the corresponding congruence for the H-sum, we apply the following identity which relates 
S- version and H- version multiple harmonic sums (see [18, Theorem 6.8]): 

l 

(-l) l &S n (s)= (-1)' (35) 

where \_t = i Si is the concatenation of Si to and the vector s = (s r , . . . , s±) is obtained from 
s = (si, . . . , s r ) by reversing its coordinates. We substitute s = ({2} b , 3, {2} a ) and n = p — 1 
in (35) and note that if I > 2, then one of Sj = ({2} c ) for some positive c and therefore 
flp-i(si) = (mod p). Thus 

(-l) a+b+1 S l p _i({2} a , 3, {2} 6 ) = -^i({2} 6 , 3, {2} a ) (mod p). (36) 

On the other hand, we have (see [19, Theorem 4.5]) 

flp_i({2}°, 3, {2} 6 ) = -H p ^{{2}\ 3, {2} a ) (mod p). (37) 

Now from (36) and (37) we get the required congruence and the proof is complete. □ 

Theorem 4.2. Let a, b be non-negative integers and a prime p > 2a + 26 + 1. Then 
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Proof. We begin by considering the case a > 1. Then from identity (26) with n = p — 1 and 
congruence (34) we have 



p-i 
F k=l 



4^ ff fc _ 1 (-2a)(-l) fc / px 

--Z^ p& (l-2piTfc-i(l)-rJ (modp). 

y k=i 

By the obvious equality 

H n (a)H n (b) = H n (a, 6) + H n (b, a) + F n (sgn(a6)(|a| + |6|)) a,kZ\ {0}, (38) 
we obtain 

S p ^({2} a , 1, {2} 6 ) = - |flp_i(2a + 26) + 4flp_i(l, 2a + 26) - ^flp_i(-2a, -26) 

+ 8# p _i(-2a, 1, -26) + &Hp_i(l, -2a, -26) 

+ 8# p _i(-2a - 1, -26) + 4# p _i(-2a, -26 - 1) (mod p). 

Applying (viii), (38), (vi) and simplifying we deduce that 

S P -i({2} a , 1, {2} 6 ) = ~H p ^(2a + 26) - ^H p ^(-2a, -26) (mod p). 

Now the required congruence follows from (i) and Lemma 3.1. The case a = is handled in 
the same way with the help of identity (25) and congruences (vi), (vii). The corresponding 
congruence for H p _i({2} a , 1, {2} b ) easily follows from relation (35). □ 

Corollary 4.1. Let p be a prime greater than 7. Then 

1,1, 4 ) = S p -!(1,2,2,2) = ^B p _ 7 (modp), (39) 
1 1 889 

S p _ 1 (l,l,l,6) = -ijJ_3 + — B p _ 9 (modp). (40) 

Proof. The first congruence in (39) follows from [19, Theorem 7.3] and the second one is a 
consequence of Theorem 4.2 with a = 0, 6 = 3. To prove congruence (40) we set a = 0, 6 = 4 
in Theorem 4.2 to obtain 

85 

S h (1,{2} 4 )e-£ h (modp). (41) 

On the other hand, using [19, Theorem 6.4] we can rewrite the length 5 sum S p -i(l, {2} 4 ) in 
terms of length 4 sums and S p -\(1, 8) as 

25 p _i(l, I 2 ! 4 ) =S P -i(l, 2, 2, 4) + 2, 4, 2) + 4, 2, 2) 

255 

+ 5 p _i(3, 2, 2, 2) - — 5 p _i(l, 8) (mod p). 

Now applying [19, Theorem 7.5] and expressing all length 4 sums in the congruence above in 
terms of the three quantities S p -i(l, 1, 1,6), S l p _i(l,8) = B p _g, and Sp_i(l, 2)5 p _i(l, 1,4) = 
— |-Bp_3 (mod p) we obtain 

1 281 

25 h (1,{2} 4 )e3S p _ 1 (1,1,1,6) + -5 p _ 1 (1,2)5 h (1,1,4) - — 8) (mod p). (42) 
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Comparing congruences (41) and (42) we conclude that 

1 1889 
5 p _i(l. 1, 1. 6) = -gVi(L 2)5 p _i(l. 1, 4) + 8) 



1_ 3 _ 1889 
54 B p~z + "gig" 

□ 



From Corollary 4.1, [19, Theorems 7.3, 7.5], and [36, Section 2] we obtain the following 
description of multiple harmonic sums of weight 7 and 9. 

Corollary 4.2. Let p be a prime greater than 7. Then 

(i) all multiple harmonic sums 5 p _i(s), i? p _i(s), of weight |s| = 7 belong to Q(2)i? p _7 
modulo p; 

(ii) all multiple harmonic sums S p -i(s), iJ p _i(s) of weight |s| = 9 can 6e written modulo p 
as linear combinations of the two quantities B p _g and B^_ 3 with coefficients in Q(3). 

In [35, Theorem 3.16], Zhao proved that 

5 p _i({l}° , 2, {l} 6 ) = 2, {l} 6 ) = (mod p) 

if a + b is even and a prime p > a + 6 + 3. In the next theorem we prove a deeper result on 
the above sums. 

Theorem 4.3. For non-negative integers a, b and a prime p > a + b + 3, 

(°a+i 2 ) B P-a-b-2 (mod p) if a + b is odd, 

2(a+6+3) 



Vi({ir,2,{i}^ 



^ (l + (-l) a ( a ^ 3 )) ( mod ^) ifa + bis even; 



( ( a t+i 2 ) B P-a~b-2 ( mod P) tfa + bis odd, 



F p _ 1 «ir,2,{i} b )^< 

, 2(a+6+3) 

Proof. From Theorem 2.2 for c = 2 we have 



pgp - a -"- 3 (l + (-l) a ( a t|2 3 )) ( mod ^) ifa + bis even. 



-Mi 1 ) >M A J J - 2^ t.a+6+2 + Z> p+1 



£a+b+2 1 
k=l k=l 



Setting n = p — 1 and observing that 

(>-') = (P - ')&» -»)-&-*) =(-!)>(! -^(D) (m „d/) 



we have 



o mi a 9 r U 6 W y^ l-pgfc(l) ^ ^-i(a + l)(l-pgfc(l)) 
VUl 1 } > 2 ' U) ) = " Z> i.a+6+2 Z> £5+1 



^.a+b+2 
k=l k=l 



= -flp_i(a + 6 + 2)+pflp_i(l,o + 6 + 2)+pflp_i(a + 6 + 3) (43) 
-flp_i(a + 1,6 + 1) +pflp_i(o + l,6 + 2) +pflp_i(a + l,l,6 + l) 
+ pflp_i(l,a + l,6+l) +pH p ^ 1 (a + 2,6+ 1) (modp 2 ). 
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If a + b is odd, then by (ii) we obtain 

5 p _i({l} a , 2, {l} 6 ) = -H p ^(a + 1,6 + 1) = (° +_^ + 2 ) B p _ a _ 6 _ 2 (mod p). 

If a + 6 is even, then from (43), (i)-(iii) after simplifying we get the required congruence. 
Considering equality (35) with n = p — 1, s = ({l} a , 2, {l} 6 ) modulo p we find 

( _ ir +6+l^_ i({1} 6 5 2j = _ flp _ l({ l } a ) 2 , {!}&) (mod p) 

and therefore 

f-5„_i({l} 6 ,2,{l} a ) (mod p) if a + 6 is odd, 
(mod p) if a + o is even. 

Similarly, from identity (35) considered modulo p 2 we obtain 

(-ir + Vi({i} b ,2,{ir) 

ee H p ^({l} a , 2, {l} 6 ) - H Hp^sjHp^) (mod p 2 ). (45) 

Sl U S2 = 

ai} a Mm 

It is clear that the product H p _i(si)H p _i(s2) has one of the forms 

ff H ({ir)I H ({ir, 2, {l} 6 ), 1 < m < a, 

or 

flp-i({l} a , 2, {1}')# p -i({1}^), < I < b. 

If a + 6 is even, then by (i) and (44), it easily follows that both products are congruent to 
zero modulo p 2 . Thus in this case from (45) we obtain 

S p -i({l} 6 , 2, {If) ee Hp^diy, 2, {1}") (mod p 2 ) 

and the proof is complete. □ 

Remark 4.3.1. It is worth mentioning that there is a notion of duality for multiple harmonic 
sums: given s = (si, S2, • • • , s r ) we define the power set P(s) to be the partial sum sequence 
(si, s± + S2, ■ ■ ■ , Si + • • • + s r -±) as a subset of {1, 2, ... , |s| — 1}. Then s* is the composition 
of weight |s| corresponding to the complement subset of P(s) in {1,2, ... , |s| — 1}, namely, 

s* = F- 1 ({l,2,...,|s|-l}\P(s)). 

It is easy to verify that (s*)* = s. Moreover, by [18, Theorem 6.7] we have that 

Sp-i(s) ee -£ p _i(s*) (mod p) 

for any prime p. In particular, for any positive integers a, b it follows that 

S p _i({2}°, 1, {2} 6 ) ee {2D 1 , 3, {2} 6 - 1 , 1) (mod p), 

Sp^({2} a , 3, {2} 6 ) ee {2} a , 1, {2} 6 , 1) (mod p), 

S p -i({l} a , 2, {l} 6 ) ee -5 p _i(a + 1, & + 1) (mod p), 

anc? fry Theorem 4-1 and Theorem 4-2 we are able to find S p -i(l, {2} a , 1, {2} 6 , 1) and S p -i(l, {2} c 
modulo a prime p in terms of Bernoulli numbers. 
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5. Congruences for "odd" multiple harmonic sums 

In this section we prove some congruences involving "odd" multiple harmonic sums. These 
results will be needed in the next sections. 

Lemma 5.1. Let s € N r . Then for any prime p > 3, 

Proof. Changing the order of summation in the definition of H n (s), with n = we obtain 

1 



H n (si,s 2 ,...,s r ) = ^2 



0< fel <...< fc ,<n^ + 1 ) S1 ---^ + 1 )" 

y I 

(-l)Sl+-+S T 



2^+-+ S r(k 1 + 1)«1 • • • (k r + l) S - 

0<k r <...<k 1 <n-l ( ' V 7 

(_l)|s| 

|s| if ra (g r , . . . , si) (modp). 



Theorem 5.1. For any positive integer m and any prime p > 2m + 1, we have 



ffp-i (2m) = 4m{2 m + l) pB P- 2m - 1 (m ° d p2) ' 



□ 



(46) 



Proof. Let n = (p — l)/2. It is easy to deduce (see [15, Lemma 2]) that 

77 n (2m) = F 2n (2m) - 

which by (i) and (iv) implies the first congruence. 

From [35, Section 2.4, Lemma 2.12] it follows that the "odd" multiple harmonic sum H n 
as well as H n satisfies the following stuffie relation: 

m 

mH n ({2} m ) = ^(-l) fc - 1 ^ ra (2fc) • H n ({2} m ~ k ). (48) 
k=i 

Now the congruence (47) easily follows from the above identity and (46) by induction on m. 
Indeed, for m = 1 it follows from (46), and if m > 1, then by (48) and induction hypothesis, 
we have 

mH n ({2} m ) = (-l) m - 1 F n (2m) (mod p 2 ), 
and the theorem is proved. □ 

Theorem 5.2. Let a, b be non-negative integers and a prime p > 2a + 2b + 3. Then 

C Cfoia o r 0l 6\ _ B p _2a-2b-3 (2a + 2b + 2\ 

S P _^({2} ,3,{2}) = 2a + i j (modp), 

fl«d({2} 8 , 3, {2}*) = ( _ 1)a+ Hi^3p« + 26 + 2\ (mQd 
2 a + l\za + ly 
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Proof. Setting n = ^rp m identity (23) and observing that 



2 

' J P\(3 P\ (2k L V 



\\ k Sl 2K2 2) V 2 2]_ = f^k 1 d n / 49 A 

we obtain 

SV±({2} a ,3,{2} b ) = -2H E - 1 (2a + 2b + 3) - AH^-i (2a + 1, 2b + 2) (mod p). 

2 22 

Now the required congruence easily follows by (iv) and (v). From Lemma 5.1 and Theorem 
5.1 it follows that 

Hp^i{{2} m ) = 2 2m H^ 1 ({2} m ) = (modp). 

2 2 

Therefore applying identity (35) with n = and s = ({2} a ,3, {2} b ) we obtain 

(-ir+ 6+1 ^ ({2} b , 3, {2} a ) = -H^ ({2} a , 3, {2} fc ) (mod p) 

2 2 

and the theorem is proved. □ 

Theorem 5.3. Let a, b be non-negative integers, not both zero, and a prime p > 2a + 2b + 1. 
Then 

S ¥ ({2}M,{2}") S — ( + jB^^! (modp), 

^w-H'i - ( -"Tr" 2) (M>'---'- < m ° d ">- 

Proof. If a > 1, then from identity (26) with n = (p — l)/2 and congruence (49) we find 

Sp=i {{2} a , 1, {2} b ) = -2Hj^i (2a + 2b + 1) - 4Hp^i (-2a, -26 - 1) (mod p). 

2 22 

Now the required result easily follows by (iv) and Lemma 3.2. If a = 0, then from (25) we 
have 

<SVl(M 2 } 6 ) = 2H E - 1 (-26- 1) (modp). (50) 

2 2 

Moreover 

(-»-!) = E WFT= E 



£,26+1 (p _ fc)2b+l 

k=l fc=(p+l)/2 ^ ' 

= E ^rr = ^-i(- 2ft - 1 )-^(- 2& - 1 ) Mp). 

fc=(p+l)/2 



Therefore, since 6 > 0, by (vi) it follows that 

1 2- 2b - 1 

g z ^i(-26-l) = -fl p -i(-2&-l)= 2b + 1 g p _ 2b _! (modp). (51) 

From (50) and (51) we conclude the truth of the first congruence for a = 0. From identity 
(35) we easily obtain 

(-lT +b S P _^({2} b , 1, {2D = #^({2}°, 1, {2} 6 ) (mod p) 

2 2 

and the proof is complete. □ 



(52) 
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Lemma 5.2. Let r, a be positive integers. Then for any prime p > r + 2, we have 
%I (r) = ^ LfJ (r) + C ~ I + k ) %i <T + k )P k 

-(- 1 ) E fc P (modp a+1 ). 

fc=o v 7 

If r is a non-negative integer, then for any prime p > 2r + 3, we have 

H g-i (-2r - 1) = 1 J +1 ff^i (2r + 1) (modp 2 ). 
Proof. Let n = (p— 1)/2, m = [p/4j . Then the first congruence follows easily from the equality 

To prove the second one, we have 

/ -i\k 



(-lfV til = v 



\2r+l 



-(2fc + ^(p-2A;) 2 

2 ^ (P - 4^) 2r+1 ~ S (p _ 2k) 2r + 1 



m 



2 ^ ( 4fc ) 2r+1 (P/( 4fe ) " 1 ) 2r+1 ( 2k ) 2r+1 (P/( 2k ) ~ 1 ) 2r+1 

= { H ^ 2r + !) + ^ m(2r + 2) ) 

+ ^T^( 2r+1 ) (modp 2 ). 
From (52) with a = 1 and r replaced by 2r + 1, we get 

2tf m (2r + 1) + {2r + 2 1)P H m {2r + 2) = (1 + 2 2? ' +1 )# n (2r + 1) (mod p 2 ), 
which after substitution into the right-hand side of the above congruence implies the result. □ 

6. Results on ^-analogues of Leshchiner's series 

In this section, we consider finite p-analogues of zeta and beta values arising from the 
truncation of Leshchiner's series (5)-(8) and defined by the finite sums 

,„„ ,„ ._3^ (-irH t -,({2}") f.^ (-ir^g t _ 1 ({2}"-') 

Cp (2™ + 2).= -^ t2(?) + 2^ ^ +2(?) , 

5 £^ (-l) k+m - l H k M2\ m ) 



m p—1 



( _ 1)l+m -,-ijj t _ i({2 ) 



m—j 
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(-l)*+"»-J( 2 fc fc )^ fc ({2} 



fc=o 

m (p-3)/2 ( _ 1) fe+m- J -^ gjfc({2}m - j: 
+ fr^ t^o 16 fc (2fc + l)2i+2 

/3f2m + n , 3 (p ^ /2 (-ir( 2 fc fc )^({2n 
^ p (2m + lj.- 4 ^ 16*(2fc + l) 

™ ^ (-l)^(^)ff fc ( { 2}"^) 

Applications of our results on congruences obtained in Sections 4 and 5 allowed us to establish 
the following congruences for the finite p-analogues of Leshchiner's series. 

Theorem 6.1. Let m be a non-negative integer and p be a prime such that p > 2m + 3. Then 

4m(l — 2~ 2m ) 

2p( p (2m + 2) = y - — '- B p „ 2m -i (mod p), 

2m + 1 

{m + l)(2m + 1) „ / i ■, 

C P (2m + 3) = 2^3 B P~ 2m ~ 3 ( mod P)- 

Note that for Cp(2), Cp(3), Cp(4), stronger congruences were proved by the authors in [26, 14]. 
For Cp(5), we recover the congruence Cp(5) = — 6/5i? p _5 (mod p) proved in [14]. 

Theorem 6.2. Let m be a non-negative integer and p be a prime such that p > 2m + 3. Then 
M2m + 1) = ; i6m(2m ; i} £p-2™-i (mod p), 
— 2m 2 + 3?ri + 2 

C p (2m + 2) = - 22m+3(2m + 3) P^ P -2 m -3 (mod p 2 ). 

Note that for /3 P (1) and Cp(2), sharper congruences were proved in [15]. 

7. Finite identities and their formal proofs 

By revisiting the combinatorial proof (due to D. Zagier) presented in Section 5 in [22], it 
is possible to obtain finite versions of identities (5)-(8). Here we derive these identities in a 
formal unified way based on applications of appropriate WZ pairs. 

Lemma 7.1. Let m be a non-negative integer. Then for every positive integer n, we have 



" (-l)*- 1 _lA (-l) n+fc+m J ff fc _ 1 ({2} m ) 

2^ l.2m+2 O 2^ 



k=l 



y^ ly 

Z»2m+3 o 



3 a (-i) m ff fc -i({2p (-ir- j g fc -i({2r^: 

l) k+m H k ^({2} m ) 



(53) 



k=i 



5f (-l) fc+m - 1 ff fc -i({2} m ) (-l) fe+m -^ 1 F fe _ 1 ({2} m -^ 
9 2^ Ui^fcx +2 Z^Z^' 



(54) 



2 ^ * 3 ( 2 *) Ut, ^ +3 ( 2 *) 
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Proof. For non-negative integers n, k, we consider a pair of functions 

FAnk) _ (-l) n+k (n-k- 1)1(1 + aUl-a) k 

^ k >~ 2(n + fc + l)! ' n ^ fc + 1 > 

CAnk) _ {-l) n+k (n-mn + l){l + a) k (l-a) k 
Gl(n ' fe) " ( (n + i)2_ a 2 )(n + A; + 1)! > 

where (c)fc is the shifted factorial defined by (c) k = c(c+ 1) ■ ■ ■ (c+ k — 1), k > 1, and (c)o = 1. 
It is easy to check that (Fi,G\) is a WZ pair, i.e., it satisfies the relation 

Fi(n + 1, fc) - Fi(n, fc) = Gi(n, fc + 1) - Gi(n, fc). 

Summing the above equality first over k and then over n, we get the following summation 
formula (see [14, Section 4]): 

N N N 

5^Gi(n-l,0) = ^(Gi(n-l,n-l) + Fi(n,n-l))-^Fi(iV,fe-l), iV e N. (55) 

n=l n=l fc=l 

Substituting (Fi,G\) in (55) we obtain 

" (-!)»-! _! * ! 3n 2 + a 2 n-l q2 

2^ n 2_ 2 2 ^n 2 ( 2n W 2 -o 2 mV 

n=l n=l V n / m=l (56) 



1 JV (_l)JV+fc . fl 2 , 

+ 2 ^ £.2^\ II ( X " ^2 J 



2 k 2 ( N )( N+k ) 11 
k=l ^ UA fc / m=l 

Expanding both sides of (56) in powers of a 2 and comparing coefficients of a 2m we derive the 
finite identity (53). Similarly, considering the WZ pair 

F 2 (n,k) = F 1 (n,k), G 2 (n, fc) = ( -^- G 1 (n,k) 

and applying the same argument as above we obtain the second identity. □ 
Lemma 7.2. For any non-negative integers m, n we have 



a (-i) fc _ I v (- 1 ) n+fc+m (fc fc )^*(^ 2 > m ) 



fc=0 



(2fc + l) 2m +! 4 ^ i 6 fc(2fc + 1) ( n + fe + 1 ) 

3 a (-ir0^({2p a a (-ir^- { 2 k k )H k {{2}^) 



(57) 



4 Z-> 16 fc (2fc + 1) 16 fe (2fc + l) 2 -? +1 

A 1 lA (-l)^flg t ({2} m ) 

(2fc + 1)2^+2 + 4 16fc(2A; + 1) 2(n+fc + +l) 

5f (-l)^(!)fl*(W - " ( -l)fc+^( 2 fc)ff fc({ 2} 
4^ 16 fe (2fc + l) 2 16 fc (2fc+ 1)^'+ 2 



m—j N 



(58) 



Proof. As in the proof of Lemma 7.1, applying summation formula (55) to the WZ pairs 
„, (—l) n+k (n — fc — 1)! /I + a\ /l-a\ 



Gi(n,k) 
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{-l) n+k {n- fe)!(2n + l) /l + a\ /1-cr 



(n + fc)!((2n + l) 2 - a 2 ) 



/I + a\ /I — a\ 

(—)»(—)»• 



and 



F 2(n, fc) = *:), G 2 (n, fc) = Gi(n, k), 



2k + 1 

we get the required identities. 



2n + l 
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□ 



8. Proofs of Theorems 6.1 and 6.2 
Proof of Theorem 6.1 Setting n = p — 1 in (53), we get 

t(2m+2)-H ( 2 m 2) ^ (-^W) 
C P (2m + 2) - H p ^(-2m - 2) - - ^ ^ 



Since 



("l) fe 



m=l 



1 1 pff fc (2) 2 



we obtain 



( p (2m + 2) = - J H p _i(-2m - 2) - i-^- # p -i({2}™, 1) 



(-1)' 



2p 

-if p _ 1 ({2} m+1 ) (modp), 

which by (i), (vi), and Theorem 4.2 implies the required congruence. 

Similarly, setting n = p — 1 in (54) and applying (59), (i) and Theorem 4.1 we get 

C P (2m + 3) = H p ^(2m + 3) - - ^ -^-^-^ 

/_ -I \m-l , i x 

= il p _ 1 (2m + 3) + (-H p _ 1 ({2} m+1 ) + flp_i({2}"\ 3)) 



\m-l 



l^T3 _ + ^ — ^- l({2} ' 3) 

(m + l)(2m + l) 



2m + 3 

Proof of Theorem 6.2 

Setting n = (p - 3)/2 in (57), we get 



B p -2m-3 (mod p). 



8<2rn + n-H i 2 m ., , (-l)<--')/2+m (-l) t fflH t ({2}") 
+ -1) + j £ 16i(2 , + 1)( ^,/ r) - 



Taking into account that 



( 2 fc fc ) 



2(2k + 1) 



fc-i 



1-2A; AAV 



3=0 



(2J + 1) 5 



2p 



2fc + 1 



(59) 



(mod p 2 ), (60) 
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we obtain by Theorem 5.1, 

fc=0 

_ f_l)(P+l)/2+m _ 
= H P _ 1 (-2m-l) + ^^- H P _^({2} m ,l) (modp 2 ), 

2 Z 2 

which implies the first congruence of the theorem by Lemmas 5.2 and 5.1, and Theorem 5.3. 
Similarly, from (58) by (60), Theorem 5.1, Lemma 5.1 and Theorem 5.2, we have 

C p (2m + 2) = Iti (2m + 2) + ^-L (h^ ({2} m+1 ) + pH^ ({2} m , 3)) 

Ptfp-2m-3-\ ~ P-TCp^i(|2} ,3) 



2 2m+3 (2m + 3) 

9 

■pBp- 2m -3 (modp 2 ) 



2 



2m 2 + 3m + 2 , , , 2 



2 2m + 3 (2m + 3) 
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